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Abstract 



We propose that the massive compact halo objects (MACHOs) observed 

in the recent microlensing experiments with an apparent best fit mass of about 

l/^ ' 0.5Mq are objects made out of "mirror" baryonic matter rather than famiUar 

OP ■ baryons. Such a possibiUty arises naturahy within the framework of mirror 

pq ■ matter models proposed recently to accomodate the sterile neutrinos that 

^^ , seem necessary to solve all the neutrino puzzles simultaneusly. We show that 

Qs I for mirror matter parameters that fit the neutrino observations, the maximum 

mass of mirror stars are of order 0.5Mq and their main sequence lifetime is 

PLh' much less than the age of the universe. They are therefore likely to be black 

Q ■ holes. Mirror matter machos have the advantage that they do not suffer 

from the problems encountered in the conventional red, brown or white dwarf 

j^ [ interpretation. The calculations also apply to the question of how the world 

of familiar matter would be different if all fundamental mass parameters were. 

X 

C^ ■ I. INTRODUCTION 

The nature of the dark matter in the universe, for which there is considerable observa- 
tional evidence, is a mystery |jl|. There are a number of experiments in progress to resolve 
this mystery. Here we will address the issues raised by the microlensing experiments ^ 
which monitor millions of stars in the neighbouring Large Magellanic Cloud (LMC) to see 
if any them suddenly brighten for a certain duration of time and then fade away again. The 
brightening observed is attributed to a lensing effect due to passage of a dark massive object, 
presumably from our galaxy, in front of the LMC star. The duration At of the brightening 
of a MACHO event has been calculated and is proportional to y/m/v where m is the mass of 
the object and v is its velocity. If one assumes that the object is from our galaxy, its velocity 
is determined; as a result, from At, one can deduce the macho mass. Based on the 14 events 
from the MACHO and the EROS collaborations that are attributable to MACHO events, 
one obtains a best fit mass [Q] of O.SMq for the MACHOs. Also it has been established 
that there are no MACHO candidates with masses between 10~'^Mq to 10~^Mq. It is also 



expected that these objects can comprise as much as 30 to 50% of the halo mass fraction 
0. The question then arises: "what are these objects?". 

The simplest possibilities are conventional baryonic objects such as red, brown or white 
dwarfs whose masses are expected to be in this range or that they are neutron stars. It 
has however been argued by Hegyi and Olive |Q that a large class of baryonic candidates 
are incompatible with observations. More recently, Freese et al ^ have made a detailed 
study of the possibility that they could be red, brown or white dwarfs and have found such 
interpretations to be highly problematic. (For a more recent analysis of some of these issues, 
see 10.) We will summarize these difficulties in a subsequent section. Accepting them for the 
moment their conclusion then leads us to search for alternative explanations. Our finding in 
this paper is that if there exists a mirror universe with identical particle and force content to 



the visible universe prior to gauge symmetry breaking |p!2| , p!0| ,p| , then for certain choice of the 
symmetry parameters, the maximum mass of the "mirror" stars is of order 0.5Mq and they 
could therefore be the machos observed in the microlensing experiments. Since they are made 
of mirror baryons, they avoid all the problems encountered by machos made of conventional 
baryonic matter [0,§ . We also estimate the main sequence lifetime of the mirror stars and 
find that in the parameter range of interest it is much less than the age of the universe. As 
a result, the machos are in the form of black holes. Let us note that precisely such models 
P] have recently been proposed to accomodate all the neutrino observations by identifying 
the lightest of the mirror neutrinos with the sterile neutrino needed in understanding the 
LSND result together with the solar and the atmospheric neutrino results. 

While this is the main result of our paper, our study also applies to the question of how 
the world of familiar matter would look if the masses are all scaled by a common factor. 
Similar questions have been studied in the past [0-|^. 

The main results of this letter are the following: we review the scaling laws for maximum 
and minimum values for the stellar masses as the masses of the electrons, W-bosons, protons 
and neutrons vary together {rrii —>■ (rrii). The maximum value, which is of particular interest 
for us is derived by using an argument in the literature that the radiation pressure inside 
a stable compact stellar object should not exceed the gas pressure inside it and by finding 
how that condition scales as the elementary particle masses vary. We find that, for ( of 
order 15, mirror stars have maximum mass of the order of the O.5M0 needed for MACHOs 
in the halo. We give qualitative arguments to support the hypothesis that the initial stellar 
mass function (IMF) for the mirror sector is likely to peak near the maximum value for ( in 
the parameter region of interest to neutrino physics. These two arguments {( ~ 15 explains 
both MACHO masses and the neutrino results) in our opinion strengthen the conjecture 
that mirror matter machos are the dark massive halo objects in our galaxy seen in the 
microlensing data. 

II. SCALING LAWS FOR STELLAR MASSES AND MAIN SEQUENCE 

LIFETIMES 

We begin with a brief overview of the mirror universe model and the the parameters 
describing fundamental forces in the mirror sector. As mentioned, one considers a duplicate 
version of the standard model with an exact mirror symmetry which is broken in the process 
of gauge symmetry breaking. All particles and parameters of the mirror sector will be 



denoted by a prime over the corresponding familiar sector symbol- e.g. mirror quarks are 
m', d' , s', etc and mirror Higgs field as H', mirror QCD scale as A' . We assume that < H' > 
/ < H >= A'/A = (. Since one expects the masses of the neutron and proton to be given 
by the scale A and charged lepton (and current quark) masses to be given hy < H >, scaling 
both parameters by the same amount implies that all fermion masses that are relevant to 
the discussion of stellar structure scale by the same amount, i.e. we have rrii —>■ C,mi with 
i = n,p, e, W, Z. Furthermore, this gives weak cross sections varying as C~^ for fixed values 
of energy. With these simple rules, assuming that electroweak and strong coupling constants 
do not change, we can say a great deal about how the properties of stars would change. 
We start with the four equations of stellar structure: 

dP/dr = -Gp{r)M{r)/r^ (1) 

dM{r)/dr = A'Kr'^p{r) (2) 

L(r)/47rr2 = -{lQ/'i)asB{T^ / pt^)dT/dr (3) 

dL/dr = A-Kr'^t{r)p{r) (4) 

where ^(r) is the opacity (cross section per unit mass) at radius r, asB the Stefan-Boltzmann 
constant, L(r) the luminosity at radius r, and e(r), the rate of energy generation per unit 
mass at radius r. We will need three terms in the equation of state (below) taken one or 
two at a time: 

P = {p/m)kT + {AasB/^c)T^ + {h'' /2m,){:i/M)\p/mfl^ (5) 

where the three terms represent gas pressure, radiation pressure, and (non-relativistic) 
degenerate electron pressure, m is the nucleon mass, irif. that of the electron. We have 
neglected such niceties as keeping track of how many objects there are for each m of gas (2 
for H, 3/4 for He, etc) 



We will make standard, illuminating if crude, approximations |14] in order to understand 



the C, behavior of the solutions to the above equations. First we write 

P = p GM/R, p = 3M/47r R^ (6) 

where P and p are roughly core averages. Here M and R are mass and radius of core or 
star; our approximations are not good enough to be precise on such points. (In practice we 
will adjust M to be the mass of the whole star and R will fall short of the core radius for 
the sun). Equation (^ gives the useful relation 

P = (47r/3)i/=^G'M2/3p4/3 (7) 

To find the minimum mass of a star, we neglect the radiation pressure term in Equation 
(^), insert into Equation (0), solve for T and maximize with respect to p, giving 

kT = (G'V2)(87r/3)2/3M^/3^«/3^e//i' (8) 



Following Phillips |T^, we set T = Tig, the lowest temperature that gives sufficient burning 



to match energy escape and solve for M, obtaining 

M™„ ~ [h^mg/{m,G'm''/')f^ (9) 

We know from more detailed analysis that Mmm is of the order of 0.07Mq and T^g ~ lO^i^. 
We use Equation @ to obtain the variation with (. m and r/ig go as (. T^g, in principle 
must be found by solving the four coupled Equations (|I]-^. Roughly, however, nuclear 
binding energies will go linearly with ( so we will approximate the variation of Tig as linear 
as well. We will see below that the solution of approximate equations gives T varying with 
( (numerically) not greatly different. We thus obtain 

M^in ~ C"' (10) 

We can also, again following Phillips fl^, use Equation (^ to find the maximum mass 



of a (main sequence) star. As the mass of the star gets bigger, the core temperature rises. 
Therefore, of the three terms in the expression for the pressure in the Equation ^ we expect 



Pg and Pr to dominate. Following Phillips |T^, we parameterize them as fractions of the 
total pressure P as below: 

Pg = PP, Pr = {l-/3)P (11) 

We eliminate T and solve Equation (|lT]) for P, obtaining 

PP = [{p k/m)\(3-' - l)/{AasBl^c)Y" (12) 

Using Equation (|^ again then gives 

M^,, ~ [(1 - (3)c/asBY"G-''\k/mf/(3' (13) 

As (3 approaches 1, the energy density is increasingly dominated by photons (relativistic 
particles) and stars become unstable. Taking a cutoff around (3 ~ 1/2 gives a maximum 
stellar mass around TOM©. Thus the range for stars is roughly O.OTM© to TOM©. From 
Equation (|1^) one sees, in the approximation that instability sets in at the same j3 indepen- 
dent of C, that Mmax varies as C~^ (like Mmin)- It is similarly easy to see from the standard 
expression for the Chandrashekar mass that it too varies as (~'^. Note that, in a model with 
r/ig varying linearly with (, but m constant, both Mch and M^ax would be independent of 
( while Mmin would go as C"^''^ (since higher mass for the electron would permit contraction 
to higher densities before Pauli repulsion becomes important). Note that such is the case for 
the mirror matter model investigated in references ^ in order to solve the neutrino puzzles. 
We now consider stellar burning as ( varies. We approximate Equation (^) as 

L = {16n/3fasB{RT)V{^M) (14) 

where k is the opacity, for which we keep just 7 — e scattering contribution i.e. we take k 
afl 



^Omitting other contributions to opacity tends to overstiniate k as C increases and hence under- 
estimate the luminosity and overestimate the main sequence hfetime. For our purpose, therefore, 
our assumption about k is a conservative one. 



K = (TT/me ~ C I^Q (15) 

Since the rate of energy generation is determined by the weak interaction rate for p+p — »• 
e"*" + (i + z/, we approximate Equation® by 

L = eM (16) 



with 

e = WvEppp/w? (17) 

We take Wv from Clayton H^ and, as above, continue in (^, obtaining 

e={pEjm')C-'f{T,/C) (18) 

where Tg = T/10^, Epp is the energy release for the full pp chain into other than neutrinos, 
and 

f{x) = 3 X 10-^^x-2/3e-3=^-^i^'^' [1 + 0.021a;-i/3 + O.Olx^/^ + 9.5 x x 10"^] (19) 

One factor of C,~^ in Equation (p!8D comes from the parenthesis preceding C~^ and two 
come from the behavior of the weak cross section, taking into account energies increasing 
with C- 

Equating the two expressions above for L permits us to solve for RT as a function of T, 
while Equations (|^) and (|^ give RT = MGm/k. Combining these results we can write, 

M' = A{3/lQiTf{k/mGy{e^Kjp^asB)T'Cf{n/0/f{15) (20) 

where we have normalized to the temperature (Tg = 15) at the center of the sun. 

Equation (|20| ) gives an approximation to the variation of stellar masses with (. With it, 
we can solve for R, p, L, and the main sequence lifetime Ims = (O.lMc^/L). In Figure 1 
(a,b) we give the results, as a function of M, for temperature T, radius R, and main sequence 
lifetime, tMs for C = 1-0, 15. We have inserted an overall factor to scale the main sequence 
lifetime of the sun (^ = 1, < Tg >~ 7.5) to 10^'' years. The scale units are lO^i^, lO^^cm, 
and 10^ years for the three quantities. It should be noted that, in the approximations made, 
the solar (core) radius comes out to about 0.3 x lO^^cm while, in real life, two thirds of the 
sun's mass (with temperature Tg > 7) extends out to about 1.5 x lO^^cm. 

We see from Figure (1) that the ranges of radii and main sequence lifetimes fall with 
increasing ( while that of temperatures increases. In Figure (2) we address this increase in 
more detail by plotting T against ( for four values of an index that runs from 1-100 as M 
varies from Mmin to Mmgx in equal logarithmic increments. Roughly, we see that T goes 
as C^^^, so that the assumption above that T^g ~ C is not grossly out of line. It should be 
noted that, for massive stars (with (^ = 1), the pp cycle on which the above considerations 
are based is replaced by CNO cycles (Clayton, ref. Dill), in which C, N, and O catalyze 



He production from H at high enough temperatures. Since the weak interactions in the 
cycle are all weak decays, which increase, in rate, linearly with (, we expect that massive 
stars, for large (, should have even shorter main sequence lifetimes than those of Figure 1. 
Thus the lifetimes in Figure 1 should be considered upper bounds. However estimating the 
abundances of C, N, O could be quite complicated (see below). 

We will use, below, the unsurprising result from Figure (1), that for C > 5 main sequence 
lifetimes are all much shorter than the age of the universe; they fall roughly as (~^. 



III. MIRROR VS BARYONIC MACHOS 

In this section we show that mirror matter MACHOs (MMMs) provide a good explana- 
tion of microlensing events within which microlensing data can determine the parameter C. 
Finally, we discuss briefly tests of the model. 

First, mirror matter resolves a number of MACHO problems. Fields, Freese and Graff [Q, 
in a very detailed work, raise several problems with baryonic MACHO candidates, including: 

-All baryonic candidates require that the MACHO population be near the minimum 
of the range permitted by observations and that Qb be near the maximum of the range 
permitted by BBNS theory if the sum of individual baryon components is to be less than 
the total baryon number. MMMs avoid this problem completely since mirror matter does 
not enter into the baryon budget. 

-Neutron stars and black holes from supernovae cannot fit the fact that lensing events 
point to MACHO masses around O.SMq. We will see below that mirror matter does provide 
an explanation of the mass observation. 

-Brown dwarf explanations conflict with a growing number of observations that show 
that the index (3 in the initial stellar mass function (IMF), N{m) ~ m~^, which is over 2 
(2.35 is a commonly accepted value |jl5|) for m > Mq whereas it is well under 2 ioi m < Mq] 
under 2 means that most of the mass is in the higher mass stars. The conflict is because 
the frequency of microlensing events appears to require a MACHO population with a total 
mass on the order of up to half Qb while the decrease in f3 for low masses precludes such a 
result. Since mirror matter is not baryonic, it has no such problem. 

-Finally, the "favorite candidate," white dwarfs, suffer from several problems raised by 
Fields, Freese and Graff 0, including: (i) not being seen - for example in the Hubble deep 
field (see Flynn, Bahcall, and Gould [|l6l) as they should be since 0.5Mq dwarfs can only 



cool slowly; (ii) a need for a large population of galactic massive stars and their supernovae 
to produce galactic winds to cleanse the galaxy of the processed material in the white dwarf 
ejecta; and (iii) a contradiction between the amount of carbon the progenitors would produce 
and the amount observed. None of these problems would arise with MMM progenitors. 

We turn now to the question of why the MMMs should have masses around O.SMq. This 
would be the case (see Figure (1)) if we have ( greater than say, 15, since the maximum 
stellar mass then falls in the region between 0.5 and 1.0 solar masses. The work of [^ shows 
in some detail that such values are just what is required to provide simultaneous solution 
of the atmospheric, solar, and LSND neutrino problems (and provide warm dark matter in 
addition). There are, furthermore, reasons why (a) the mirror matter stellar IMF would 
peak near the maximum mass and (b) remnant masses would be similar to initial masses. 
Both stem from the decrease in cross sections with increasing (. 



Current theories of star formation (see, for example, Adams and Fatuzzo |]T5| and ref- 
erences therein) from molecular cloud core collapse require a mechanism to stop accretion 
during collapse, i.e. to limit the size of the star. Such mechanisms are based on scattering, 
but cross sections for scattering of molecules, atoms, ions, and electrons off photons, atoms 
and molecules will fall as (~'^. Thus it is not unreasonable to expect that the mirror matter 
IMF should become more and more strongly peaked near Mmgx as ( increases. Additionally, 
we might expect a modest increase in ("^ M^ax in the Section 2 estimate since scattering pro- 
cesses that create instability as /9 ^ will become less effective leading to smaller /? values 



and, through Equation ([T3|), larger values of Mmax- 

Cross section decrease also predicts that there should be little mass loss between the 
initial star and the remnant. As ( increases, the radius of the star decreases and neutrino 
cross sections decrease. Thus neutrino confinement times decrease sharply and it is doubtful 
that (Type II) supernova shocks can be formed, thereby creating the population of MMax 
black holes that are detected as MACHOs. These two C > 15 features - decreased mass loss 
and mirror star masses peaking around 0.5Mq - result in fixing ^ in a range optimal from 
the point of view of fitting neutrino results in reference [§]. 

In sum, MMMs appear to have a number of positive features as the explanation of 
microlensing events. 

Finally, we turn to the question of observational tests of the MMM hypothesis. The 
following come to mind: 

-Absence of any optical observations of lensing objects as lensing events accumulate; 

-Within the qualitative picture above, relatively strong peaking of lens masses into a 
narrow range; 

-Possible detection of the black holes by some new method; unfortunately estimates by 



Heckler and Kolb [0 show that, even with new instruments (the Sloan digital sky survey 
telescope), black holes under IOMq could saturate the halo mass without being detectable 
from the signal from interstellar material infall. 



-Possible future detection of black hole MACH binaries ||T8|] that if they exist in sufficient 



number through the emission of gravitational waves in experiments such as LIGO, VIRGO, 
TAMA and GEO. 



IV. LUCKY TO BE ALIVE 

Finally we note some of the implications for the familiar world from the above results 
on varying (. As noted in Section I, there is a growing literature on the changes that 
would obtain if standard model, and other parameters, were different 10^ • The present 
investigation adds to those results. 

The two most important changes as ( grows would appear to be the absence of supernovae 
and the decrease in stellar lifetimes. These imply that, as ( grows, there would be lower 
abundances of heavy elements in the interstellar medium with which to form planets and 
carbon based life forms on them, as well as less time in orbit around main sequence stars 
during which it would be possible for the latter to occur. Although rates for radiative 
processes would increase linearly with (, main sequence lifetimes would fall as (~^ as shown 
in Figure 1. It is the factor of p^ in phase space that, apparently results in increasing rates 
and decreasing cross sections. 

As ( decreases, the combination of decreasing rates and increasing cross sections would 
be likely to interfere with current models of star formation cited above. For example, collapse 
times ([Gp]~^/^) would increase like C~^ while cross sections for scattering that would tend 
to disperse the collapsing cloud would increase like C~^- 

In conclusion, we have studied the variation of stellar mass and lifetimes as the masses of 
the elementary particles rrie, rrip, rrin, Mw all vary in the same way (given by the parameter 
Q. We conclude that for a value oi C, — 15, the maximum mass of the mirror stars is around 
half a solar mass; as a result, they could be viable candidates for the MACHOs observed in 



the various inicrolensing searches^. The many problems encountered in trying to explain the 
0.5Mq machos white dwarfs and brown dwarfs etc are now easily avoided. We do a crude 
analysis of the dependence of the "main sequence" life time of the mirror stars as a function 
of the ( variable and find that for similar ( values, the mirror macho is most likely a black 
hole since main sequence stellar lifetimes are few times shorter than the age of the universe. 
We also propose several tests of this hypothesis. 

We are grateful to D. Clayton, S. Nussinov and D. Rosenbaum for many useful dis- 
cussions. We are grateful to K. Freese and W. K. Rose for reading the manuscript and 
comments. The work of R. N. M. is supported by the National Science Foundation grant 
under no. PHY-9802551 and the work of V. L. T. is supported by the DOE under grant no. 
DE-FG03-95ER40908. 



^Note that, in the mirror matter model of references |^,^ the neutron is unstable for (" in the 
range of interest, so stellar structures become either mirror white dwarfs or black holes. Our 
considerations with respect to the maximum mass of such objects would still obtain. The result 
would be prompt black hole formation and the application to the MACHO problem would be the 
same as in this letter. 
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Figure Caption 

Figure 1 (a, b): Temperature T, radius R and the main sequence lifetime tMs of the mirror 
stars as a function of the stellar mass M for six different values oi ( = 1.0, 15. The units for 
the above are lO^K for T, 10^° cm for R and 10^ yrs. for tMs- 

Figure 2. Variation of temperature T as a function of ( for 4 of 100 equal mass steps between 
Mmin (step 1) and Mmax (step 100). Step 3 for ^ = 1 corresponds to the sun. 
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